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ABSTRACT 

We describe the classical cosmological tests, such as the LogiV-LogS', redshift- 
magnitude and angular diameter tests, and propose some new tests of the 
evolution of galaxies and the universe. Most analyses of these tests treat the 
problem in terms of a luminosity function and its evolution. The main thrust of 
this paper is to show that this is inadequate and can lead to a incorrect conclusions 
when dealing with high redshift sources. We develop a proper treatment in three 
parts. 

In the first part we describe these tests based on the isophotal values of the 
quantities such as flux, size or surface brightness. We show the shortcomings of 
the simple point source approximation based solely on the luminosity function 
and consideration of the flux limit. We emphasize the multivariate nature of 
the problem and quantify the effects of other selection biases due to the surface 
brightness and angular size limitations. In these considerations the surface 
brightness profile, and the distribution of the basic parameters describing it, play 
a critical role in modeling of the problem. In general, in the isophotal scheme 
the data analysis and its comparison with the model predictions is complicated. 
In the second part we show that considerable simplification is achieved if these 
test are carried out in some sort of metric scheme, for example that suggested 
by Petrosian (1976). This scheme, however, is limited to well resolved sources. 
Finally, we describe the new tests and compare them to the traditional tests 
demonstrating the observational and modeling ease that they provide. These new 
procedures, which can use the data to a fuller extent than the isophotal or metric 
based tests, amount to simply counting the pixels or adding their intensities as a 
function of the surface brightness of all galaxies instead of dealing with surface 
brightnesses, sizes and fluxes (or magnitudes) of individual galaxies. We also show 
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that the comparison of the data with the theoretical models of the distributions 
and evolution of galaxies has the simplicity of the metric test and utilizes the data 
as fully as the isophotal test. 

Subject headings: cosmology:theory-galaxies:evolution-galaxies: luminosity 
function-galaxies:photometry. 
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1. INTRODUCTION 



Galaxies and other extragalactic sources provide the most direct means of studying 
evolution in the universe. This is done using the classic cosmological tests such as the 
angular diameter-redshift-magnitude relations or the source counts (also known as the 
LogiV-LogS 1 ) test. (For a general description of these tests see, e.g. Weinberg 1972.) These 
tests, which rely primarily on the distribution of the magnitudes or fluxes of the sources, 
have had limited success in determining the cosmological parameters and/or the evolution 
of galaxies. There are two fundamental reasons for this failure. The first is the well-known 
difficulty of disentangling the evolution of the sources (see e.g. Tinsley 1968 or Tinsley 
& Gunn 1976) from the evolution of the universe (Weinberg 1972). As a result, over the 
years, the focus of such studies has been shifted from the determination of the cosmological 
parameters to the evolution of galaxies in different assumed cosmological models. The 
second difficulty arises from the fact that galaxies are extended (i.e. resolved) sources 
and there is ambiguity in defining proper magnitudes (or luminosities L and fluxes /) 
and diameters. In addition the samples of sources are not merely limited by their fluxes 
(or magnitudes) but there exist other selection biases or data truncations due to surface 
brightness or size limitations. These aspects of the problem are usually ignored. This may 
be an acceptable approximation for high surface brightness sources at low redshifts, but it 
is woefully inadequate when dealing with data at high redshifts extending to low surface 
brightness sources. The main purpose of this paper is to present a proper analysis of the 
various observational biases that are encountered in this process. There are two ways one 
can carry out this task. From the original data sample one can select a subsample with 
fewer and simpler biases (as we do in §3), or one may correct the model expectation fully for 
all known selection biases, which is the approach we take in the rest of the paper. The first 
method is more appropriate for tests dealing with the moments of the distributions such as 
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flux-redshift or size-redshift relations. The second method is preferable when dealing with 
the various source counts and uses all of the valuable data. 

The bias due to the magnitude or flux limit is accounted for by various means. The 
most common practice is to use isophotal values, i.e. the values of these quantities up to or 
at some limiting apparent surface brightness &i S0 . However, because of the rapid decline of 
the apparent surface brightness b (defined as flux per unit angular area of a resolved source) 
with redshift z (see, e.g. Tolman 1934), 



where B is the intensity or the absolute surface brightness at the source, the biases due to 
the surface brightness and size limits of the observations become important at high redshifts 
and/or for low surface brightness sources. These effects are often ignored or are dealt with 
indirectly by using a limited portion of the available data. 

The corrections required for these effects, sometimes referred to as aperture corrections, 
inevitably require the knowledge of the surface brightness profile 



and the distributions of the central surface brightness B , the characteristic or core radius 
r , and other parameters ctj defining the profile f(x). For example, a = lor 1/4 for disks or 
spheroids, respectively, where the profile is described by the simple relation ln/(x) = — x a . 
Early examples of methods to correct the redshift magnitude relation for the aperture 
effect were described by Sandage (1972), using an iterative procedure, and by Gunn & Oak 
(1975), assuming a fiducial cosmological model. These authors used empirical relations for 
the luminosity within the radius r, 



b = B/(l + z)\ 



(1) 



B{r) = B f(r/r ) with /(0) = 1, 



(2) 
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It was shown by Petrosian (1976) (P76 hereafter) that these corrections can be carried out 
more directly. It was also shown in P76 that one can separate the evolution of the surface 
brightness Bq from the evolution of universe, and can avoid some of the above difficulties by 
dealing with the angular sizes and magnitudes up to and within a "proper metric" radius 
r p obtained from a specified value of the quantity 

Fix) n dhar . , r , iS 

r l = -^77^ = 2 m t( \ Wlth x = - 4 
x z j{x) dmL{r) r 

which is equal to the ratio of the average surface brightness within r to the surface brightness 

at r. 

The above equations describe the source brightness profile and its basic properties in 
terms of two convenient parameters; the central surface brightness Bq and core radius tq. 
These parameters are not easily accessible to observations and their relative values for 
different values of a are difficult to interpret. This difficulty can be overcome if we transfer 
the above relations to observationally more meaningful parameters. One commonly used 
such sets of parameters is the effective radius and surface brightness. The effective values 
refer to the radius containing half the total luminosity L which means F(r e g/ro) = 0.5F(oo). 
The ratios of the effective to central values are: B eS /B = 0.189,2.54 x 1CT 2 ,3.45 x 1CT 3 
and 4.66 x 1(T 4 , and r efr /r = 1.66, 13.5, 1.82 x 10 2 and 3.46 x 10 3 for a = 1, 1/2, 1/3 and 
1/4, respectively. 

Figure |1] shows the profile /, the curves of growth of luminosity F/F(oc) and the 
function 77 as a function of (r/r e g) a , for the above values of a. This shows the general and 
relative characteristics of these functions. 

To demonstrate the effects of the redshift, in Figure ^| we show the variation with 
redshift of the fraction of the luminosity and the area (expressed in magnitude units) 
within a surface brightness limit b for a = 1 and 1/4 profiles. Instead of the surface 
brightnesses B and b we use the more familiar magnitudes; /1 = — 2.51og6 + const., and 



Fig. 1. — Variation of surface brightness profile /, normalized luminosity growth curve 
F/F(oo), and the inverse of the i] function defined in equation (HI) versus {r/r e ^) a for a = 
1, 1/2, 1/3, and 1/4 (solid, dotted, short dashed, and long dashed lines, respectively). Note 
that In/ = — (r /r cS ) a (r cS /r ) a and the effective radius r eff is defined as F(r/r cS ) = 0.5F(oo). 

/i e fj = — 2.51og(£> c ff) + const.. As evident the observable area and luminosity vary rapidly 
with redshift, specially for low surface brightnesses and have a different behavior for the 
two profiles. This will produce a variation with redshift of the relative abundances of disks 
and spheroids. Another possible representation of the above graphs will be in terms of the 
proper metric radius r p defined above. This is preferable because this definition of radius 
relies on the data within some measured isophot and not on the unobserved outer parts 
which are needed to determine the effective radius r cS . This procedure will be developed 
further in §3. 

In this paper we review several old procedures and propose some new ones for the study 
of the evolution of galaxies, and possibly that of the universe, whereby instead of dealing 
with individual galaxies we deal with the combined brightness of all galaxies. The new 
methods simplify the data analysis enormously and are perfectly suited for modern digitized 
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data. In §2 we first give a brief description of the proper analysis of the classical tests 
for isophotal quantities that includes all the observational selection effects as well as the 
effects of the surface brightness profile, and treats the problem in terms of the multivariate 
distribution ip(B , ro, ctj, z) instead of the commonly used luminosity function cf)(L,z). 
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Fig. 2. — Variation of the fraction of the luminosity (solid lines) and projected surface area 
(dashed lines), in some arbitrary magnitude units, within a surface brightness isophot /i 
versus redshift for several values of the effective surface brightness // e ff- a) for an exponential 
(disk) profile, In/ = — x; b)for a de Vaucouleurs profile, In/ = — x 1 ^ 4 . Note that for each 
value of n e Q the luminosity and area shrink rapidly as we approach the maximum redshift 
Zmax = (Bo/b) 1 / 4 — 1, where B and b are related to fi c g and /i as described in the text. 
Also note that this decline is more pronounced for disks than spheroids indicating that the 
relative populations of sources with different profiles will vary strongly with redshift. 

There are several reasons for the popularity of the latter procedure. The first is that 
we have accumulated a considerable knowledge about the local luminosity function of 
galaxies (see e.g. Efstathiou et al. 1988 and Lin et al. 1966) but little information on 
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the distributions of the B and r (see, e.g. Sandage & Perelmuter, 1991). Secondly, 
until recently it was believed that that the distribution of the surface brightness was 
fairly narrow (Fish, 1964; Freeman, 1970). However, in recent years, because of increasing 
interest in low surface brightness dwarf galaxies, some data has been accumulated on the 
intrinsic distributions of these parameters indicating broad distributions (see e.g. McGaugh 
1996; de Jong 1996; Tully & Verheijen 1997; Impey & Bothun 1997 and references cited 
there). Of course, one can use the luminosity function by replacing the total luminosity 
L = (47r)7rro-B F(oo) [see eq. (|3])] for either B or r . In any case a multivariate description 
is required. In §3 we repeat the analysis of §2 for metric quantities. In §4 we describe the 
new tests and their relations to the multivariate distribution i[) and the profile f(r). Finally, 
in §5 we give a brief summary. 



The classical test use the observed relations between the magnitude (or flux /), angular 
size ( radius 9 or area irO 2 ) and redshift to determine the cosmological parameters and the 
evolution of sources as described by the general "luminosity" function tp. The cosmology is 
introduced via the relations 



where di and 6a = g?l/(1 + z) 2 are the luminosity and the angular diameter distances, and 
Qi represent the cosmological parameters such as the density parameter Q, the deceleration 
parameter go o r the cosmological constant A (see, e.g. Weinberg 1972). All the classical tests 
can be described in terms of the observed distribution of flux, size and redshift; n(l,9,z). 
For the purpose of the demonstration of the effects that we would like to emphasize here, 
let us consider the cumulative source counts as a function of redshift, which we denote by 



2. 



CLASSICAL TESTS: ISOPHOTAL VALUES 




(5) 
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N(> /, > 9,z). The cumulative and differential counts of the so called LogA^-LogS 1 relation 
is obtained by the integration of the above expressions over the redshift. 

For simplicity, in the above relation and in what follows, we ignore cosmological 
attenuation, if any, assume either bolometric (or monochromatic) fluxes, so that we can 
ignore the K-correction, and assume spherical symmetry with brightness profile f(r) 
independent of the wavelength. The complications due to K-correction, asphericity, etc. can 
be easily included in the relations that follow. We will address some aspects of these in §5. 



2.1. Point Sources 



The usual practice (see, e.g. Metcalfe et al. 1995 or Tyson, 1988) is to compare the 
observed cumulative and differential distributions iV(> l,z) and n(l,z) = — <9iV(> l,z)/dl 
(and their integrals over redshift) with that expected from models via the relations 

dV r°° dV 
N(> l,z) = — : <f)(L,z)dL and n(l, z) = — 4ird 2 L (l)(4ird 2 L l, z), (6) 

dz JAnd'j^l dz 

where V(z, O^) is the co-moving volume up to redshift z and p{z) = J °° <ftdL is the co-moving 
density of all sources at z. Here and in what follows we assume a complete 47r sterradian 
sky coverage. These expressions is what one would expect for unresolved or point sources, 
where only the flux limit counts in the selection process (see, however, a modification below 
in §2.3). 



2.2. Extended Sources 



For extended sources such as galaxies the selection process is more complex and 
additional corrections are required. We now describe these selection biases. 



2.2.1. Surface Brightness Limit 



A source to be detected must have an apparent central surface brightness exceeding 
the detection threshold which must be several times the standard deviation a of the 
fluctuations in the background brightness. We denote this limit by &i S0 . If we ignore the 
image degradation due to the finite size of the instrumental and atmospheric point spread 
function (PSF), which can be done if the core size r 3> 9 s d,A, and if the pixel size is 
less than or comparable to the width 6 S of the PSF, then the surface brightness selection 
criterion requires that 

B > B iso = (1 + z)\ so . (7) 

However, for small sources or high redshifts the effect of the finite size of the PSF cannot 
be ignored and the selection bias is more severe than indicated by this relation. For a 
PSF= g{6/6 s ) the surface brightness is modified to 

B{r) = B f(r/r , 8 s d A /r ), where f = f * g (8) 

is the convolution of the the actual profile with the PSF. As a result, the central surface 
brightness is reduced by £ = B /Bq, where 

Z{0,d A /r O ) = J f(8d A /r o )g(6/e s )6d0/ J g(6/6 s )6dd. (9) 

For the purpose of illustration let us consider a box PSF with the radial width of 9 S . This 
reduction factor then simplifies and equation (|7|) is modified to read 

B > B iso x 2 s /F(x s ) with x s = 9 s d A /r . (10) 

Similar expressions can be derived for other forms of the PSF. For i s < 1 this reduces to 
equation (^) but its effects become important for x s near unity i.e. for partially resolved 
and unresolved sources. Note that 6 S is replaced by the pixel size # P i x if 9 S < 6 pix 



2.2.2. Size Limit 



Another criterion for selection of extended sources such as galaxies is that their sizes 
must exceed some limit. One way to quantify this is to have the isophotal angular radius 
(namely, the radius where the surface brightness has dropped to the specified isophotal 
value &i SO ) be larger than some specified size 9. If 9 ^> 9 S , then the isophotal angular radius 
is given by 

9 iso = (r /d A )f-\B iso /B ), (11) 

where the function f" 1 is the inverse of the profile function /. Then the selection condition 
$iso > 9 is satisfied if 

B > B iso /f(9d A /r ). (12) 

However, as 9 decreases toward 9 S one should use the modified profile / of equation (g) 
in place of /. In any case, it is clear that for 9 > 9 S this inequality will provide a more 
restrictive limit than the surface brightness limit because it requires that more than one 
pixel to exceed the surface brightness limit b\ so . This can be demonstrated mathematically 
also by setting 9 = 9 S in the last equation and comparing it with the limit in equation 
(|T0|). The ratio of the two limits is equal to rj(x s ) which according to equation (§) is greater 
than one, except for the unlikely event of the surface brightness increasing with r. In the 
opposite case when 9 < 9 S one is dealing with unresolved or point like sources in which case 
a size limit does not make sense. 

For non-spherical sources we can follow a similar procedure by dealing with the isophotal 
angular area uo (which for spherical sources is equal to ir9 2 ) as the area of the sources with 
apparent surface brightness b > b lso . However, the relation of this to the surface brightness 
profile will be more complicated. For example, for elliptical sources with a constant 
ellipticity we can express the profile / as a function of the area a/a^ with a = 7rr 1 r 2 , 
where r\ and r 2 are the core radii along the major and minor axes. For randomly oriented 
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elliptical sources this will amount to replacement of the quantity flc^/Vo in equation ([12]) by 
(ujd 2 A / do) 1 ! 2 . The distribution function ip now will be a function of a and the ellipticity or 
the ratio r\jr2- 



2.2.3. Flux Limit 



Finally the sample of sources is subject to a flux limit. In this section we consider the 
flux Z iso within the isophotal angular radius or up to the surface brightness limit 6 iso . 
The flux limit then implies that l- lso = L{9 mo dA) / ^d 2 L is greater than some specified flux I. 
Using equations @ and (||) we can write this limit as 

B F(x iso ) > d 2 L l/{ixrl) with x iso = 6 iso d A /r = f-\B iso / B ). (13) 

Note that the left hand side of this inequality is independent of r . However, if 9 iso ^> 9 S , 
then one should replace the profiles / and its integral the luminosity growth curve F by 
the corresponding values, / and F, modified by the PSF. In this case the above relation 
becomes more complex with the involvement of the additional variable 9 s dA/ro (see §2.4 
below). 



2.2.4. Combined Limits 



Thus, for a given value of 9, /, 6; so and 9 S the above three inequalities determine the 
region of the B — r plane that is accessible at these particular conditions. This region varies 
with redshift becoming smaller at higher redshifts. The redshift dependences are hidden in 
.Biso and d&- Figure |3] show the three boundary conditions obtained by the equality sign 
in equations flTO), (jT2|) and flT3| ) for the exponential, a = 1, and de Vaucouleurs, a = 1/4, 
profiles. We plot B /B iso versus r /(9 s dA), which is valid at all redshifts. The lowest 
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(heavy) solid line shows the truncation due to the surface brightness limit. The lighter solid 
lines show the effects of the size limit for several values of 9/9 s . The dashed lines show the 
truncation due to the flux limit for different values of the ratio I / (it 9 2 b iso ) . Sources lying in 
the region above all three lines are the ones which satisfy all the selection criteria. It is clear 
that, as long as 9 > 9 S , or I > n9 2 b iso , which obviously will be the case for resolved sources, 
the surface brightness limit due to the PSF, described by equation ([H]), is never important. 
However, both size and flux limits could be important depending on the relative values of 
the observational limits. For larger values of the ratio l/(n9 2 bi SO ) the flux limit provides 
the major constraint. In the opposite case the size limit becomes more important, and as 
evident from the above figures for I = Tt9 2 b- lso only the size limit is relevant. 



2.2.5. Extended Source Counts 



We can now relate the observable N(> I, > 9, z) to the distribution function ip(B , r , z) 
by the integration of the latter over the accessible region as determined by the observational 
limits. For the general case this gives 

dV ( f B 0,cr roo roo roo \ 

N{> l,> 9,z) = — / dB dr ^{B , r , z) + dB rfr ^(5 , r , z) . 

aZ \JB iBO Jro,i JBo, cr Jrofi J 

(14) 

Here r 0) i and r ,2 are obtained by solving equations (|13D and (|T2| ) for r in terms of B and 
other observables, and 

B , cr = B iso /f (rT 1 (l/(n9\ so ))) , (15) 



is the intersection point of the the two boundary conditions described by equations (|T 
and (p~3[ ) or the intersection of a solid and a dashed line in Figure |3|; r/ _1 is the inverse 



function of the function rj defined in equation (|j). Note that as stated above for I = Ti9 2 b- ls 
this critical value of surface brightness becomes equal to -Bi SO , the first double integral in 
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Fig. 3. — Central surface brightness B in units of -Bi SO = h so (l + z) 4 versus the size r 
in units of 6 s d A depicting the surface brightness limit (the heavy or lowest solid line); 
the angular size limit for 9/9 s = 1,2,5 and 10 (thin solid lines); and the flux limit for 
l/n6 2 s b mo = 1,2,5,10,100 and 1000 (dashed lines). a)for an exponential (disk) profile, 
In/ = — x; b)for a de Vaucouleurs profile, In/ = —x 1 ^. Note that the solid and dashed 
lines crosss each other where the limiting factor changes from size to flux. The region to the 
right of the appropriate sold line before these intersections and to the right of the appropriate 
dashed line after the intersections are accessible to observation. Note also that, because B iso 
and (1a increase with redshifts, the observable region shrinks systematically, moving towards 
the upper right hand corner, as redshift increases. 
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the left hand side of equation (|H|) vanishes and we are left with the size limited part of this 
expression only. 

It is clear, therefore, that the relation for the counts of extended sources is considerably 
more complicated than the commonly used relation (|6|). In order to see these differences 
more clearly we can rewrite the above expressions in terms of the luminosity L. For example, 
if we replace r by L = (A^ttTqBqF '(oo) we can rewrite all the boundary conditions in 
terms of B and L instead of B and r . The three limits in equations (fLOf ), (0) and (|T3|) 
now give, respectively, the conditions 

B >B iso (F- 1 (X)) 2 /X, with A = 4vr4(vr^ iso )F(oo)/L, (16) 

L > L min , e = 47r4(vr^6 iso )F(oo)/ (xf so /(x iso )) , (17) 

and 

L > L min j ee 4vr4/F(oo)/F(x iso ). (18) 

Here F' 1 is the inverse function of the luminosity curve of growth F and a;i SO is a function of 
Bq/B- 1so (eq. [0). Figure f| shows the truncations produced by the above selection criteria 
in the Bq — L plane. 

If we define the distribution ip(B ,L,z) = ip(B , L, z)dr /dL, equation (fill) then 
becomes 

dV ( r B o,cr roo /-OO /-OO _ \ 

N(>l,>6,z) = — [ dB dLij(B ,L,z) + dB dLip(B ,L,z)\ . 

min,0 / 

(19) 

So far we have expressed our results in terms of the central surface brightness B Q and 
core radius r . As mentioned in §1 these parameters are not convenient for compring the 
results with observations. This task can be carried out more readily if we express the above 
relations in terms of observationally more meaningful parameters such as the effective 
radius and surface brightness defined in §1. In Figure |5] we show the size and flux limits 
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Fig. 4. — Same as Fig. |3] except for the B — L plane. The curves are plotted for the same 
limiting values as in Fig. [5]. The luminosity unit L iso = Ar[d 2 L {T[6 2 a b\ so ). 

in the surface brightness-luminosity plane for both profiles, where instead of the surface 
brightnesses &i SO and B we use the more familiar magnitudes; /ii SO = — 2.5Log&i SO + const., 
and fJ, e s,z — — 2.5Log-B c ff + 10Log(l + z)+ const.. This figure demonstrates that the region 
of B — L plane accessible to observations is different for the two profiles and shrinks with 
increasing redshift. 



2.2.6. Comparison With The Point Source Approximation 



There are several ways that this correct description differs from the approximate 
expression given by equation (fj]). Some of these were discussed by Yoshii (1993). The first 
difference is the existence of the second set of integrals in equations flTjj ) and (|B^) which 
we discussed above. Even at high values of the flux limit, I ^> 7T0 2 bi SO , when this additional 
term is negligible there are two other important differences. The first is due to the presence 
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Fig. 5. — Same as Fig. § except we plot the combined size and flux limits in terms of 
the effective values of the surface brightness in magnitudes; /z iso = — 2.5Log(6 iso )+ const, 
and [M e fi,z — — 2.5Log(-B e fj) + 10Log(l + z)+ const.. For clarity, only curves for 6 = 6 S and 
I = n6 2 s bi SO are shown. 
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of the ratio F(oo)/F(xi SO ) in the lower limit of the luminosity L m[n i, which is absent from 
the lower limit in equation The second effect is due to the breadth of the distribution 
of Bq. For narrower distributions this effect is smaller and disappears for a delta function 
distribution of Bq. 

These differences can be seen in Figures 4a and 4b as follows. For a given value of I 
or the ratio l/n9 2 b iso the point source approximation given by equation (Q) truncates the 
Bq — L plane by a vertical line at the asympthote of dashed line appropriate for this ratio. 
Then counts all the sources to the right of this line (L > 47td 2 L l) , irrespective of their surface 
brightness B , and uses <p(L,z) = f °° ip(B , L, z)dB . Equation (|19"D, on the other hand, 
indicates that this is an overestimation of the counts and that we should count only sources 
which lie to the right of both a solid and a dashed lines appropriate for the ratios 9/9 s and 
1 1 '7r0g&is O , respectively. It is, therefore, clear that ignoring the surface brightness and size 
limitations can cause significant errors the extent of which can be quantified only if we know 
the distribution function ■?/>. We have a good knowledge of the dependence of -0 on L at low 
z and high Bq but we have only scanty information on its form at high z and low values of 
Bq. The primary aim of the cosmo logical tests under the discussion here is to determine the 
variation with redshift of the general distribution function ip. A detailed investigation of 
these aspects are beyond the scope of this paper. Here we make some simple comparisons 
between the extended and point source results which do not require a knowledge of the 
distribution of Bq. 

1) Luminosity limits. The ratio of the limiting luminosities L min j and/or L m[n to the 
limit 47rd|Z of the point source approximation, which depends only on the surface profile, 
are shown in Figure for a = 1 and 1/4, respectively, and for several values of the surface 
brightness (actually the ratio Bq/B 1so ) and the ratio g = I / (n6 2 b- lso ) . We use the effective 
rather the central surface brightness and express the above ratios in magnitudes. As evident 
for high values of the ratio g, i.e. for higher flux limits, there is an increasing bias against 
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detection of extended sources at higher redshifts. At lower values of this ratio there is 
additional bias against detection of galaxies at low redshifts due to the size limit. 
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Fig. 6. — The ratio of the minimum detectable luminosity for extended sources in the 
isophotal scheme [eqs. ( fP7D and (|l"8Dl to that of point sources versus redshift, for indicated 
values of the effective surface brightnesses (more exactly the ratio Bo/b lso ) and for four 
values of the ratio g = l/{j0 2 bi SO ) = 1,3.16,10 and oo, represented by long dashed, short 
dashed, dotted and solid lines, respectively. a)for an exponential (disk) profile, In/ = — x\ b) 
for a de Vaucouleurs profile In/ = — x l / A . Note the rapid suppression of the flux, as in Fig. 
2, when we approach the maximum redshift z max = (i?o/&iso) 1//4 — 1, where the relation of Bq 
and &i SO with // e g and /i; so , respectively, are described in the text. Note also the differences 
between disks and spheroids as in Fig. |[ 

2) Redshift Distributions. These differences can also be seen when we compare the 
redshift distributions expected for point and extended sources. For the purpose of 
illustration let us assume that L and Bo are uncorrelated (clearly not a good assumption) 
so that we can separate the distribution function as -0 = p(z)h(B )<f)(L/L*)/L*, where p(z) 
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and L*(z) describe the density and luminosity evolution of the sources, and h(B ) gives 
the distributions of the central surface brightness. If we define the cumulative functions 
H(B ) = / °° h(B )dB and $(ar) = /~ <p(x')dx', with H(0) = $(0) = 1, then equations © 
and flT9l) become, respectively 

N(>l,z) = ^-p(zM^d 2 L l/L*) (20) 

and 



N(> I, >9,z) = —p(z) f °' C " h{B )${L min ,i/L*)dB + r h{B )${L min> o/L*)dB ) . 

dz \JB mo JB , cr J 

(21) 

In absence of an exact knowledge of the distribution h(B ), we compare these expressions 

for different assumed values of the central or effective surface brightness which amounts to 

a delta function approximation of h(Bo). In this case equation ( pTj) simplifies to 

dV { $(L minj/ /L*), if B 0iCr > B ; 
N(>l,>e,z) = —p(z)\ (22) 

Assuming a Schechter luminosity function, <p{x) oc x p e~ x , we evaluate the redshift 
distributions for some representative values of the surface brightness Bo (or p e s) and for 
several combinations of the limits I (or magnitude m), 6 and 6i SO (or p\ so ). The results are 
shown in Figure [7] for p — — 1, and for a cosmological model with Q — 1 and A = 0. We also 
assume absence of any density or luminosity evolutions; p and L* constants. 

It should be noted that the relative shapes of the point source and various extended 
source distributions are independent of the cosmological model or the density evolution p(z). 
As evident equation @ or (|2CiD give quite incorrect redshift dependences, overestimating the 
number of sources by a large factor at high redshifts, specially for low values of the surface 
brightness due to the surface brightness limit, and at high values of surface brightness due 
to the size limit. Clearly ignoring these effects could lead to incorrect results. For example, 
if these expressions were used to derive the extent of the luminosity evolution, L*(z), they 





Fig. 7. — The redshift distribution N(> I, z) of equation^]) of the point sources (top 
heavy solid lines) and extended sources, equation (^), for seven different values of the 
effective surface brightness. All curves are normalized by the peak value of the point source 
distribution. In each case the dotted lines are for low surface brightnesses (n — — 1, 1 and 
3) falling below the turning point of the solid lines in Fig. 4c, where the surface brightness 
limit is most important, and the dashed lines are for high surface brightnesses (n = 7, 9 and 
11) above this turning point, where the size limit is important. The thin solid line (with 
n = 5) lies generally near the turning point. The limiting magnitude is m = 2.51ogZ + const 
and m iso = — 2.51og(7r6^6 iso ) + const.. The surface brightness is in units of magnitude per 
square arcsec. a and d for disks, and b and c for spheroids with different limits. 
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would underestimate this evolution by factors equal to F(x iso )/F(oo) and x 2 so f(x iso )/F(oo), 
depending whether l/(n9 2 b iso ) ^ 1 or is equal to 1, respectively. The situation is more 
complicated when the effects of the dispersion of the surface brightness or its correlation 
with the luminosity (or core size vq) are included. In such cases there would be errors in the 
determination of the density evolution as well. 

3) Source Counts. Integrating the redshift distributions over z gives the cumulative 
counts. Differentiation of this gives the differential counts. For example, for point sources 
n(l) = Jo°° ^(j)(4:7id 2 L l, z)(A r Kd 2 L )dz. Figure |] shows the differential magnitude counts of 
extended sources with various values of surface brightness as well as that of point sources. 
We use the same model parameters as above. Again, as evident the neglect of the selection 
aspect discussed above can cause considerable error in the determination of the evolution 
of the general luminosity function or the cosmological parameters. 



2.3. Point Sources Revisited 



It should be noted that some of these effects are present even for unresolved or point 
sources. However, in this case the correct equation is only slightly different than equation 
(|) or (|2"0l) . When 6 s dA S> ^efr the sources are unresolved and they all have essentially the 
same profile g and size 9 S as the PSF; B{9) = Bo£ t (9 s dA/i r o)g{9/9 s ). The total luminosity 
can be written as L = And\{n6 2 s )B^. If we limit ourselves to the isophotal fluxes > / 
and sizes > 9, then the three selection criteria (eqs. [fTO], [I2| and |13fl) become almost 
identical in form; B > B- mo / (Ki^(9 s dA/ro)), where £ is defined in equation (|9|) and 
Ki = l,g(9/9 s ) and (47r)(7r^6 iso )G'(6 l i SO /6' s )//, respectively for the three limits. In the above 
relations we have defined the cumulative PSF as G(x) = Jq 2g(x)xdx, G(oo) = 1, and #i SO is 
obtained from g{6\ so /6 s ) = B- lso /(Bo£) = L iso /L. Clearly the size limit does not make sense 
for unresolved sources and the flux limit is the most restrictive limit. It can be shown then 
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Fig. 8. — The differential counts of point sources (top thick line) and that of extended sources 
with different surface brightness limits (thin lines with n — 1 to 6 from top to bottom) for 
the same model parameters as Fig. |7|. a) for an exponential (disk) profile, In/ = —x; b) 
for a de Vaucouleurs profile In/ = — x 1//4 . Note the reduction at large magnitudes for high 
surface brightness sources and large suppression at all magnitudes for low surface brightness 
sources. 
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that the correct expression for the source counts is 

A\T poo poo 

N(> l,z) = — dB / dLiP(B , L, z), (23) 

where the lower limit of the luminosity is obtained from the solution of 

LG (g-\L- lso /L)) = Andll, with L iso = 4Trd 2 L (n6% so ). (24) 

Here g~ x is the inverse function of the PSF. In addition to the surface brightness cutoff 
and the integration over Bo, this expression differs from the simple equation ([]) also by 
the presence of the term involving the cumulative PSF G in the integration limit. This 
difference becomes important only for flux limits very near the isophotal values; I = n6 2 b\ so . 
The surface brightness limit can be important for unresolved galaxies because of their low 
intrinsic surface brightness or effective temperatures. But for other point sources such 
as quasars whose surface brightness is equal to that of a hot accretion disk this effect is 
negligible (becoming important only at extremely high redshifts) and the point source 
approximation of equation @ is very accurate. 

2.4. Combined Counts 

In principle we can combine the counts of the resolved and unresolved sources by 
replacement of the B(r) with the modified profile B(r) of equation @. We can then repeat 
the procedure carried out for the extended sources with the replacement of the profiles / 
and F with / and F and change the limits correspondingly, except now the profiles are 
functions of the additional parameter OgdA/ro. However, now the size limit is unnecessary 
because we can include all sources. Of course one must make sure that the sample of 
galaxies, for example, is not contaminated by other unresolved sources (e. g. stars). We 
therefore have the simpler expression 

A\J roo poo 

N(> l,z) = — dB I dL^(B 0l L, z), (25) 

CtZ J B lso J L m i n ,l 
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where L min j is given by equation ( p"8D with -F(x iso ) — > F(x iso ); note that F(oo) = F(oo). The 
parameter x iso now is obtained from /(x iso , 6 s d,A_/r ) = B iso /B with txt\ = Lj (4irB Q F(oo)) . 

2.5. Other Tests 

Differentiation of N(> I, > 8, z) gives the differential distribution n(l, 8, z) from which 
we can calculate various moments and compare them to observations. For example, the 
flux-redshift relation can be obtained from 



< 



l(z) >= J™ (LF(x iso )/ (4tt4F(oo))) n(l, > 8, z)dl j N{> /, > 8, z). (26) 



In a similar fashion one can derive <8>—zot<8> — <1> relations. 

3. CLASSICAL TESTS: METRIC VALUES 

Some of the complications evident in the above analysis can be avoided if instead of the 
isophotal sizes and fluxes we deal with some metric values of these quantities. For example, 
if we define a proper metric size r p corresponding to a constant value of the function rj, 
say r/o, as defined in P76 and equation (||), then the expressions for the surface brightness, 
size and flux limits (or truncations) become considerably simplified. For a limiting surface 
brightness fei SO , a limiting angular radius 8 and a limiting flux I these truncations are 
described by the following: 

Bo > B iso /f(Q, r > d A 8/( } and L > 4^ d\l F '(oo) / F((), (27) 

where r](() = i] and £ = r p /r . 

Note that in contrast to the complicated truncations we found for the isophotal case, 
equations fllTf ), (12) and fllTf), the current truncations are much simpler; they depend only 



-27- 



on one observational limit and the redshift. The above limits are good for 6 6 S , the size 
of the PSF. This will always be true in this case because of the need to have a well defined 
surface brightness distribution. 

If the truncation limits are chosen so that / = 7r# 2 &i SO , then the last limit in equation 
(f2~7j) due to the flux limitation falls below the other two and can be ignored. In this case 
the data truncation in the B — r plane is parallel to the axis making the calculation of 
the observables straight forward and free of the complex limits of integration. Since the 
flux limit I does not enter in the determination of the observed distribution of the sources, 
such a sample will not be appropriate data for tests based on source counts as a function 
of the flux I. But such a sample can be used to obtain the distributions of the angular size, 
average surface brightness or redshift. For example the latter is simply 

(PJ roo roo 

N((,z) = — dB dr ^(B Q ,r ,z). (28) 

dz JB lso /f(0 Jed A /c 

Because of this simplification, such samples are well suited for tests based on the moments 
of the observed distributions. For example the angular size-redshift relation is simply given 

as 

C dV c°° c°° 

< n92 ^ >= at,s ~T~ \ dB ° / dr Trrfy(B , r , z). (29) 

N(Q, z)d\ dz JB iso /f(0 Jed A /c 

Similarly, for the flux-redshift relation we have 



< l ^ >= Mfs U2 T- / dB ° B » / dr Trr&(B , r , z), (30) 

N((, z)d l L dz JB lBO /f(0 Jed A /( 

or, in terms of the luminosity L 

<l(z) >= F ^/ F ^ ^ T dB o r dLL$(B ,L,z), (31) 

N(C, z)d 2 L dz JB lBO /f(Q JL min>( 

with a similar expression for the size- z relation. The integration limit 

£ mi n,c = 47r4(7r^6 iso )(F(oo)/C 2 )( J Bo/Aso). (32) 
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These expressions are considerably simpler than the corresponding equations for the 
isophotal analysis. 

Similar expressions can be derived for other definitions of the metric quantities. For 
example, as mentioned in §1 in connection with Figures 1 and 2, instead of r p one can use 
the effective radius r e g, within which resides a certain fraction (usually half) of the total 
light. This would amount to a new definition of the constant ( as F(Q = F(oo)/2. This 
may be a convenient procedure for nearby galaxies but not at high reshifts, because it 
relies on the knowledge of the total flux, the determination of which lies at the heart of the 
difficulty associated with these tests. The procedure proposed in P76 relies only on the data 
within a specified radius and not on the data from the outer, invisible parts; 8 P is obtained 
by setting the ratio of the average to limiting surface brightnesses to a fixed value. 

4. NEW TESTS 

The discussions in the above two sections demonstrate that an accurate analysis of the 
extragalactic data for the purpose of cosmological tests is complicated and must include all 
of the above mentioned considerations. In particular, it is imperative to keep in mind the 
multivariate nature of the problem and to account for the surface brightness limitation, (eq. 
0), common in all of the above expressions. The dispersions in the distributions of B , r 
or L, and the correlations between these can have substantial effect on the final results. 
These effects are more pronounced when dealing with the isophotal quantities than with 
the metric ones. However, the latter procedure must be limited to well resolved sources, 
while the former, in principle, could be extended to unresolved sources if a good knowledge 
of the PSF is at hand. 

This task, however complicated, can be carried out given the knowledge of the 
distribution function and the brightness profile. With sufficient care in the analysis of the 
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data and in modeling one can determine either the cosmological evolution of the sources 
(i.e. the redshift variation of the distribution tp) or the cosmological parameters. Such 
analyses, which may be simple or appropriate for data limited to low reshifts is not the 
simplest method to determine the cosmological or galactic evolutions. For example, the 
traditional method of identifying sources with some apparent flux may not be necessary or 
be the most straight forward way of carrying out this task. The complexities described in 
the previous sections are the result of the multiple selection criteria needed for counting 
individual galaxies. 

We now describe two new and much simpler tests which combine the good aspects and 
avoid the complexities of the two methods described above, and are much better suited 
for the analysis of modern digitized CCD data. The essence of these tests is to reduce the 
selection criteria to one, namely the surface brightness, and deal with the distribution of 
Bq. In practice this amounts to simply counting the number of pixels at a given surface 
brightness (or adding up their intensities) independent of which galaxy they belong to. 
This way one can avoid the complexities arising from the need to define the sizes and fluxes 
(isophotal or metric) for every galaxy in the field. 

4.1. Sky Covered By Galaxies 

The first of these tests, which is related to the angular diameter test, involves 
computation of the fraction of the sky which is covered by all galaxies above (cumulative) 
or within (differential) a given range of surface brightness, b to b + db. Observational 
determination of this fraction is simple. It is accomplished by counting the number of pixels 
with a given intensity value. The expressions relating this quantity to the cosmological 
models and to the properties of the galaxies are decidedly (somewhat) simpler than those 
for the isophotal (metric) treatment. Let us first consider well resolved sources, namely 
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those with angular radii > 6 which we take to be ^> 9 S . In this case the only data truncation 
arises from the size limit which is same as that described by equations (ITT]) and (12) with 
the isophotal quantities replaced by those for an arbitrary value of b. 

B z =(l + z)\ x b = e b d A /r = f- 1 (B z /B ) with B > B z f f(9d A /r ). (33) 

Now following the same steps as in the previous sections the sky fraction covered by 
all galaxies down to a given apparent surface brightness b is obtained by adding the 
contribution i\Q\ = Tr(x b r /dA) 2 of each galaxy. The result is 

^ sky (> b, z) = — / dB x 2 b dr Tirlip(B ,r , z), (34) 

CLZ J B z JOcIa/xi, 

where the line element 

dx/dz = (dV/dz)(47Td A )-\ (35) 
Alternatively, in terms of the luminosity distribution 

fly roo -t;2 

F sky (>b,z) = -± dB — b 
dz Jb z B Jr 

where 



.F sky (> b,z) = -± dB / dLI4{B Q , L, z), (36) 

dz Jb z B q F{oo) JL min!b 



L min , b = 47r4(7r^6)F(oo)/ (x 2 J(x b )) . (37) 

Note that these equations have the simplicity of the tests based on metric values; equation 
(|34T), except for the term x\ is identical to (< 7i9 2 (z) >)N((,z), shown in equation (|29|). 
However, more importantly, the data analysis is enormously simpler because it does not 
require determination of the surface brightness profile and the metric or isophotal values for 
each galaxy. 

As in the case of the classical tests, this test also can be carried out for unresolved 
sources. But this is not much different than counting sources because all unresolved sources 
have essentially the same area, tt$ 2 . Following the procedure in §2.3 it can be shown that 
the sole truncation due to the surface brightness limit, B £ > B z , is equivalnt to 

L = And 2 A {ne 2 s B )i > L b = A^d 2 L {n6 2 s b) (38) 
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and the actual angular radius of each source is given as b /9 s = g 1 (B Z /B Q ^) = g 1 (L b /L). 
Thus, the fraction of the sky covered by unresolved sources then becomes 

F sky (>b,z) = —(7r6 2 s )j B dB j L dL(g-\L h /L)) ^(B ,L,z). (39) 

Similarly, as in §2.4, if we have a good knowledge of the form of the PSF we can 
combine resolved and unresolved sources by expressing the above relations in terms of the 
modified profile / as the convolution of the profile / and the PSF g [eq. (0)]. Then we do 
not need not to specify a size limit and the only truncation comes from the value of the 
apparent surface brightness b. However, in this case we have the added complication due to 
the dependence of the characteristics on the ratio 9 s dA/ro. For example, equation (|36|) now 
becomes 

dy roo dBn r°° - 
f skY (> b,z) = -± - / dLLx 2 b ^(B Q , r , z), (40) 

dZ JB Z B F(00) J L^b 

where now L mirijb is given in equation (|37|) with / — > / ', F — > F and x b = 9 b dA/r is 
a function of both Bq and L (or tq) and is obtained from the inversion of the relation 
f{x b ,e s d A /r ) = B z /B £, with 7rr 2 = L/{B F{oo)). 

The differential distribution can be obtained from f skj (b,z) = — <9jF sky (> b,z)/db. The 
integration of either distribution over z gives the differential or cumulative distributions of 
all galaxies irrespective of their redshift. 

4.2. Total Sky Brightness 

The second test, which is related to the flux- redshift test, deals with the contribution of 
all galaxies to the sky brightness within a range of (or above) a given surface brightness b. 
This amounts to adding all of the intensity values of the appropriate pixels. This is to be 
then compared with the expression for the total intensity (flux per sterradian) as the sum 
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of the contribution l{d b ) = irrQB F(xt>)/ (47rd|) of all resolved galaxies with 8 b > 8. 

T sky (> 6, z) = 7~r~~ ~J~ r dB »lP\ r dLL ^( B o, L, z). (41) 
(1 + zp dz Jb z -t (oo) JL min , 6 

Note again the similarity of this expression to equation ([31]). As above, if we redefine the 

profile as the convolved / and include all resolved and unresolved galaxies in the analysis, 

we obtain the relation: 

W> 0, = J^f z 12 *B a dLL^m, L, z). (42) 

The differential distribution is obtained as i s ^ y (b,z) = — <9Z sky (> b,z)/db, and the integrals 
of these over z give the cumulative and differential total sky brightness as a function of b. 
Note that i s k y (b, z) = bf^ib^z). 



4.3. Average Sky Brightness 



The ratio of the quantities described in above two tests gives the average surface 
brightness at or down to some surface brightness b. This quantity as expected is 
independent of the cosmological model parameters and depends only on the surface 
brightness profile and redshift and consequently, as already pointed out in P76, can be used 
to to determine the evolution of the surface brightness. In general, this relation is more 
complicated than envisioned in P76 where the discussion was aimed at the brightest cluster 
galaxies. For a larger and varied sample of sources this relation is more complex and not as 
obvious. If, for purpose of illustration, we assume that Bq and tq or L are not correlated, 
then from equation ([10) and its counterpart for X, this ratio becomes 

< b(z,b) >=b( F(x b )h(B )dB / / x 2 b f(x b )h(B )dB ) , (43) 

\JB Z JB z J 

where H{Bq) describes the distribution of the central surface brightness. Note that for 
a delta function, or a relatively narrow distribution, the above expression simplifies to 
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< b > /b = r](x b ) as is the case with individual galaxies. This demonstrates that the surface 
brightness profile, or the function rj based on it, plays a central role in cosmological studies 
of extended sources. 

With minor modification the above expressions will be valid for elliptical sources with 
constant ellipticity. For such sources with major and minor core radii n and r 2 , the area 
within any isophotal limit is proportional to 7rr 1 r 2 , so that if we define r = (r 1 r 2 ) 1 / 2 the 
above expressions would apply, but we now have the additional integration over the possible 
dispersion of the ellipticities or the ratio r 1 /r 2 . 



5. SUMMARY 



In this paper we deal with the analysis of the distribution of redshift-size-flux (or 
magnitude) data on extrgalactic sources, in particular galaxies, which is often used for 
testing the evolution of sources and/or the universe. The usual practice is to describe 
the characteristics and evolution of the galaxies in terms of a simple luminosity function 
(f)(L, z), as if the galaxies are point sources and the data is simply flux or magnitude limited. 
We emphasize that in reality these tests are more complicated. The proper analysis must 
involve the variation of the surface brightness profile and the multidimensional distribution 
function of the parameters that describe this profile, such as core radius, central surface 
brightness, luminosity, etc.; ip(B , r , L, a, z). Neglecting these copmlexities and the 
truncatios of the data produced by the surface brightness and angular size limits can lead 
to grossly misleading results. 

There are different ways one can account for these effects. The most efficient use of 
the observations comes from comparison of the full set of the data with model predictions 
which include the effects of all biases that are encountered in the observational selection 
processes. Use of the isophotal values down to lowest posssible isophot and size is a good 
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example of this approach. We have described the correct analysis of such a data in terms of 
the multivariate luminosity function tp. We derive the relevant expressions, in terms of the 
surface brightness profile of spherical sources, which is to be compared with observations of 
the isophotal values of the fluxes, sizes and surface brightnesses. In general, because more 
than just the usual flux (or magnitude) limit enters in the analysis these expressions are 
relatively complex. Truncation of the data due to other selection effects such as angular 
extent and surface brightness thresholds come in and can play the dominant role in defining 
the content of a sample of sources. The surface brightness profile of the galaxies plays a 
pivotal role in these calculations. 

A second approach would be to select a subset of the data which yields to a more 
straight forward comparison with models. For example, if we select the more limited 
sample of large and well resolved sources, we can use fluxes and sizes related to a metric 
(instead of isophotal) radius, such as the radius defined in P76 for a constant value of the 
77 function. We derive the relevant expressions for this case and show that they are much 
simpler than the isophotal ones, and resemble more closely the simple expressions for point 
sources. This method, therefore, would be more appropriate for evaluation of the moments 
of the observed distributions in tests like the redshift-flux or angular size-redshift tests. 
The isophotal method is more appropriate approach for tests based on source counts and, 
in principle, can be used for samples of sources which include resolved as well as unresolved 
sources. With the knowledge of the PSF at hand one can use a modified surface brightness 
profile as convolution of the actual profile and the PSF. 

The above tests, aside from the complications in the modeling, suffer from the additional 
shortcoming of needing elaborate procedures for the analysis of the data. To overcome 
some of the difficulties in both of these areas we propose a new method for the analysis 
and modeling of extended extragalactic sourses for the purpose of determining either their 
evolution or the cosmological parameters. This method is very well suited to modern 
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digitized data, and amounts to counting the numbers of pixels or summing their intensity 
values. It is capable of using all the data as in the isophotal case but the expressions 
relating the data to models are considerably simpler and are similar to the metric ones. 
Their simplicity stems from the fact that they deal with the surface brightness limit alone 
and do not include selection based on fluxes or sizes. As already shown in P76, tests based 
on surface brightness tend to be more robust and simpler. 

To reiterate, these new methods clearly have the following advantages. 

a) The data analysis is obviously considerably simpler. 

b) The expressions relating the observations to the distributions of the basic properties 
of galaxies is evidently simpler: Compare, e.g. equations fll4|) and (fSj). 

c) The dependence on the cosmological parameters is also considerably more straight 
forward: Instead of the dependence on the volume V(z), and on the luminosity and angular 
diameter distances, and d^, the new tests depend primarily on the redshift and the much 
simpler line element dx/dz. For example, for models with zero cosmological constant this is 
equal to (c/H )(l + z){Vtz + 1) -1 ^ 2 , where Vt is the density parameter and H is the Hubble 
constant. 

d) Because we are dealing with surface brightness, to first order the resulting expressions 
are independent of the weak gravitational lensing effects due to the inhomogeneities of the 
intervening matter distribution (dumpiness due to galaxies and clusters). 

As mentioned above the modern digitized CCD data are ideally suited for the task 
proposed here. However, several conditions are required for the proper application of all the 
tests proposed here. A good knowledge of the background sky brightness due to all other 
sources except extragalactic sources underconsideration is needed because we wish to go to 
as low a surface brightness as possible. Application of the new methods to integrated (over 
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the redshift) values of jF sky and T sky to the whole data say, in a CCD frame, will require 
accurate flat fielding. However, for redshift dependent type analysis the requirements are 
similar to galaxy based analysis, because redshifts are known for galaxies as a whole and not 
for indvidual pixels. For this kind of studies some identification of the pixels with galaxies 
is necessary so that we can assign redshifts to pixels. Secondly, when combining resolved 
and unresolved galaxies it is important that contaminations due to other unresolved sources 
such as stars is kept to an acceptable value. And Finally, in all these tests one must deal 
with several profiles; disks, spheroids and possibly a continuum of superposition of both. In 
future works we hope to use these new methods on the data from the Hubble Deep Field 
(see Williams et al., 1996) as well as those from ground based observations. 

This work was carried out while I spent a three month sabbatical leave at the Space 
Telescope Science Institute. I would like to thank the director Robert Williams and the 
Staff of the Institute for their support and hospitality. I also would like to thank Drs. H. 
Ferguson, M. Vogeley and M. Fall for valuable discussions on the general topic of this paper. 
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